Effective "Penetration Depth" in the Vortex State of a d-wave Superconductor 

M. H. S. Amin 1 * M. Franz 2 and Ian Affleck 3 ' 4 
1 Department of Physics, Simon Fraser University, Burnaby, BC, V5A 1S6, Canada, 
2 Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA, 
3 Department of Physics and Astronomy and Canadian Institute for Advanced Research, 
University of British Columbia, Vancouver, BC, V6T 1Z1, Canada 
4 Institute for Theoretical Physics, University of California, Santa Barbara, CA93106-4-030, USA 

(February 1, 2008) 



The temperature and field dependence of the effective mag- 
netic penetration depth (A e ff) in the vortex state of a d-wave 
superconductor, as measured by muon spin rotation (/iSR) 
experiments, is calculated using a nonlocal London model. 
We show that at temperatures below T* tx y/B, the linear 
T-dependence of A^ crosses over to a T 3 -dependence. This 
could provide an explanation for the low temperature flatten- 
ing of the curve observed in a recent pSR, experiment. 

Recent experiments on quasiparticle response in the 
vortex state of the cuprate superconductors have indi- 
cated quite unexpected behavior. Krishana et al. dis- 
covered an anomalous plateau in the longitudinal ther- 
mal conductivity k xx of Bi 2 Sr2CaCu208+<5 (BSCCO) 0] 
and YBa 2 Cu 3 7 _5 (YBCO) [§ at high magnetic fields 
which they attributed to the opening of a second gap of 
d xy symmetry. Scanning tunneling spectroscopy (STS) 
on YBCO compounds || lent support to this scenario 
by suggesting the existence of localized quasiparticles 
in the vortex core regions which is possible only for a 
gapped excitation spectrum Subsequent measure- 
ments of Aubin et al. j^] discovered hysteretic behavior 
in the thermal conductivity which is a signature of the 
influence of the vortex lattice and impurities on n xx . Fur- 
thermore, at temperatures below IK they detected an in- 
crease of k xx with the magnetic field instead of decrease 
||. This is actually in agreement with the existence of 
nodes in the superconducting gap and in contradiction 
with d x 2_ y 2 +id xy symmetry scenario. Theoretically, the 
possibility of a field induced d xy gap was discussed by 
Laughlin and others 0, and attempts have been made 
to explain the observed strange behavior of k xx without 
invoking mixed order parameter symmetry 

The question of the existence of the second order pa- 
rameter was recently raised once again by Sonier et al. in 
their /zSR experiment ||. At high magnetic fields, they 
observed a flattening of A^ (defined in these experiments 
as the width of magnetic field distribution) at low tem- 
peratures in contrast to the T-linear behavior expected in 
a d-wave superconductor ■ If one assumes that A^ 

is proportional to the superfluid density p s , then such a 
flattening could be indicative of opening of a gap in the 
quasiparticle excitation spectrum, which would result in 
exponentially activated behavior of p s (T). On the other 
hand, experiment finds strong dependence of \ c s(T — ► 0) 
on the magnetic field. This argues against the gap since 



within conventional BCS type models opening of a gap 
should not affect the value of p s {T — > 0). In this letter we 



argue that the simple relation A, 



cir 



oc p t 



is not valid for 



the penetration depth extracted in these experiments at 
finite fields and that with proper definition the observed 
behavior of A^ (T) can be explained very naturally by a 
nonlocal London model for a d-wave superconductor in 
which the superfluid density p s remains linear in temper- 
ature. 

High T c materials are extremely type II superconduc- 
tors in the sense that their coherence length £ is much 
smaller than their penetration depth A. Thus naively 
one would expect that the local London model well de- 
scribes their magnetic behavior. However, in order to 
study the electromagnetic response of a d-wave supercon- 
ductor it is necessary to define a momentum dependent 
coherence length, £p = vf/ttA-p, which diverges along 
the node directions. This divergence gives rise to nonlo- 
cal dependence between the supercurrent and the vector 
potential. Kosztin and Leggett |Q showed that in the 
Meissner state, this nonlocal relation can produce a In- 
dependence of the penetration depth, instead of linear 
T-dependence, below some crossover temperature given 
by Tj£ L = Ao(£o/Ao) where Ao is the maximum gap, 
Co = vf/ttAq and Ao is the London penetration depth. 
Such an effect is field independent and can only occur at 
very low temperatures. Therefore, it cannot directly ex- 
plain the /iSR observation assuming that A e ff is the same 
as the Meissner penetration depth. We will show below 
that an analogous calculation with a new definition for 
the penetration depth which is similar to its definition 
in /iSR experiments gives Agg ~ T 3 for T < T* . Here 
T* = A (Co /d) , and d cx ^JWjB is the distance between 
vortices, with <I>o being the flux quantum. T* is now field 
dependent and could be much larger than T^ L . 

the effect of nonlocality and 



We have studied |13 



nonlinearity due to the field induced excitations at the 
gap nodes |lrj| ], in the vortex state of a d-wave super- 
conductor using a generalized London model. In Ref. 
|]l5f we defined an effective penetration depth A ff so as 
to closely correspond to the quantity measured in ^SR 
experiments. Fig. [I] presents the result of nonlinear- 
nonlocal calculation of A c ff based on the theory devel- 
oped in Refs. |p^ , p^| . The effect of the nonlinear cor- 
rections on A c ff, as shown in Fig. [I], saturates at high 
fields and stays effectively field independent for B > IT. 
Most of the field dependence of A G ff, especially at high 
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FIG. 1. Comparison between our zero T nonlinear-nonlocal 
London theory (with Ao = 1078 A) and the recent fj,SR data 
M extrapolated to T = 0. 



fields, therefore comes from the nonlocal contributions. 
The effect of the nonlinear corrections is just to shift the 
value of A c ff by a constant, which can be compensated by 
rescaling Ao- In Fig. ^, we compare our theory with the 
/iSR data reported in Ref. § . We find that it is possible 
to fit the experimental data to both nonlinear-nonlocal 
and nonlocal-only curves, with fairly good agreement, by 
just changing the scale. At high fields both curves pro- 
vide good fits. Even below IT the agreement between the 
experimental data and the nonlocal-only curve is fairly 
good; although including nonlinear corrections enhances 
the agreement especially at 0.1T. The excellent agree- 
ment between our theory and the experimental data sug- 
gests that the same effects might be important for the 
temperature dependence of A c ff , and might be in fact re- 
sponsible for the fiSK observation. 

The effective penetration depth in the vortex state has 
been also studied numerically by Wang and MacDonald 
p7{ within a lattice model of a d-wave superconductor. 
This approach is capable of correctly accounting for the 
vortex core physics, which is beyond the scope of our 
London model. On the flip side, because of the system 
size limitations, the approach of Ref. Jl7| is limited to rel- 
atively high fields B > 10T. The present approach is well 
suited to address the low-field regime (0.1T< B < 10T) 
which is of greatest experimental interest in cuprates. 

We now present a brief discussion of analytical and nu- 
merical calculations leading to the temperature and field 
dependence of A c ff. In order to avoid unnecessary com- 
plications we only explicitly consider the nonlocal effects. 
Nonlinear corrections are included only by rescaling the 
value of Ao , as discussed above. A more detailed descrip- 
tion of the theory is given in Refs. ||3Jl|. To linear 
order, the relation between the supercurrent j and the 
vector potential A in a superconductor can be written in 
Fourier space as 




FIG. 2. Temperature dependence of A cff 2 /A 2 at different 
magnetic fields. Here Ao = A(T = 0). 



where Q(k) is the electromagnetic response tensor. Ap- 
plying the linear response treatment of Gorkov equations 
generalized for an anisotropic gap, one can calculate the 
kernel Q(k) in weak coupling limit. One finds jl4| 



Q«(k) 



4ttT 




(2) 



where 7k = Vi?-k/2, Ao = \] c? I A-ixe^VpN-p is the London 
penetration depth, lu 71 = nT(2n — 1) are the Matsubara 
frequencies and the angular bracket means Fermi surface 
averaging. Eq. (Q) is valid for an arbitrary Fermi surface 
and gap function. 

Besides the temperature dependence contained in the 
Matsubara frequencies in (||) , the gap itself has temper- 
ature dependence which becomes important at tempera- 
tures approaching T c . Assuming an isotropic Fermi sur- 
face, we use the simple form Ap(T) — A(T) cos 29, for the 
gap where 9 is the angle between the internal momentum 
p and the x-direction. In weak coupling limit and in the 
absence of magnetic field, the function A(T) can be ob- 
tained by solving the ordinary BCS gap equation which 
in d-wave case is 



v 



cos- 
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ej + A(T) 2 cos 2 29' 



(3) 



j k = -(c/47r)Q(k)A k . 



(1) 



where V is the interaction potential in the d-wave chan- 
nel and e p is the dispersion relation. At finite magnetic 
fields this equation will change and a full self consistent 
calculation will be necessary to obtain A(T). However, 
since the magnetic fields of our interest are far below the 
upper critical field H C 2 , and since we focus only on the re- 
gions outside the vortex cores, we can assume (to a good 
approximation) that Eq. (||) holds even in the presence 
of a weak magnetic field. 

One can write the London equation purely in terms 
of magnetic field, by eliminating j from (|l|) and using 
Ampere's law j = (c/47r)V x B 
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kx [Q-^kJCkx B k )] =0. 



(4) 



In the vortex state it is furthermore necessary to insert 
a source term F(k) on the right-hand side of the Eq. (^) 
which accounts for the phase winding around the vortex 

cores. We use a source function with the usual Gaussian 
cutoff F(k) =e"«o fc2 /2 m_ 

As in Ref. JIB], we define A e ff by 



SB 2 



(5) 



where SB 2 is the second mom ent of the field distribution 
in the vortex lattice and SBq is the same quantity for the 
magnetic field Bo(r) obtained by solving the ordinary 
London model on a triangular lattice with the same B 
and Ao- We emphasize that this way of defining A c ff is 
roughly equivalent to the way it is computed from the 
^SR data. Using (Q) and (||) we get 
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cif 
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k _^ (1 + £ijhkj) 2 
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where k are the reciprocal lattice wave vectors and 



A,(k) 



det Q(k) 



c*- 1 



x e 



^ fc 4 

k#0 



(6) 



(7) 



Notice that only wave- vectors of 0(d^ 1 ) and larger con- 
tribute to the vortex lattice field distribution and the 
response to a magnetic field at these wave- vectors can be 
much different than at zero wave-vector in a supercon- 
ductor with gap nodes. At low T, only the regions near 
the gap nodes are important in the calculation of Qij (k) . 
We can therefore linearize the gap as Ap w Ag = 2A(T)9 
where 9 is the angle measured from the node, and let v> 
in Eq. (Q) only take the node directions. As a result, Qij 
will be diagonal in the 45° rotated frame. 

Defining qi.2 = (k x ± k y )/y/2 we can write Qij = 
\o 2 [l + K(qi,T)]Sij, We then have 



£ n (<l)=X 2 [l-K(q 2 ,T)} 
£2 2 (q)=Ag[l-^(gi,T)] 
£i 2 (q) =£ 2 i(q) =0, 



(8) 



This is actually a Taylor expansion in K which is small 
at low B and T. Substituting into (Js[) we get 



A c ff 
A^ 



C ^ K(q u T)q 2 + K(q 2 ,T)q 2 



q#0 



where q = qd, C = Cd A and K{x,T) = K{x/d,T). We 
have also taken the upper cutoff Q" — d/£o to infinity 
since it does not affect our calculations (we keep £<j finite 
in our numerical calculation but our results are insensi- 
tive to its exact value). 
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FIG. 3. A comparison between our nonlocal calculation 
and jiSK data. 



At T = for qC, d < 1 one can write @ K(q, 0) 
-(■jt 2 /8)Q q- Substituting back into (j^) we get 



rJAeff(T = 0) OC Cd 



(10) 



This is in complete agreement with our numerical calcu- 
lation and /iSR data j| (c.f. Fig. |l|), and also with other 
calculations flij . 

At T > 0, Eq. (|) can be written as 



SX 



off 



A 



C ^ 5KjquT)g + 5K(q 2 ,T)g 2 (n) 



q#0 



g 6 



where SX cS = X(B,T) - X(B,0) and 6K(q,T) = 
K{q, T) - K(q, 0). We follow Ref. |1 by writing 



SK(q,T) = K(0,T)F 



(12) 



where t = T/T* , K(0,T) w -2(ln2) T/A is the kernel 
obtained in local approximation, and F(z) is a universal 
function which can be approximated by 



F(z) Kl-CiZ 
F(z) » c /z 2 



for 
for 



z < 2 
z > 2. 



(13) 



Here, Co and c\ are constants. Substituting ( p"5| ) and 
(|l|) into @ we find that for t < 1, falls into z > 2 
regime for all reciprocal lattice vectors q. This immedi- 
ately gives 5X c ff oc T 3 . At higher temperatures on the 
other hand, there will exist significant number of points 
with z < 2 which would give linear T behavior. Thus in 
general 



SXaS oc T 3 
<5 Acs ex T 



for T < T* 
for f«T« T c . 



(14) 



Fig. U shows the results of our numerical calculation of 
A e ff as a function of temperature. We take Aq = 1078A, 



3 



which is the value that produces the best fit to the exper- 
imental data (cf. Fig. 3), and n — A /£o = 68. We also 
assume a triangular vortex lattice aligned with x and y 
(a and b) directions. We have tried changing the core 
size as well as the shape of the vortex lattice and our re- 
sults remained unchanged. The penetration depth is no 
longer linear at low temperatures; unlike the superfluid 
density (the upper curve) . The deviation from linearity is 
stronger at higher fields in complete agreement with the 
/iSR observation || and previous theoretical work [p"7| . 
At high enough fields, the T dependence is in complete 
agreement with the T 3 -form obtained analytically. The 
curves for different fields join to a single curve at higher 
temperatures. This is essentially because the nonlocal 
corrections are most pronounced for the quasiparticles 
close to the gap nodes. At higher temperatures the re- 
sponse becomes dominated by the quasiparticles far from 
the nodes which feel much weaker nonlocal effects. 

Fig. H presents a comparison between our results and 
the experimental data. In order to get a good fit, we 
set A = A(0) = 2.65T C which is also what one obtains 
from the gap equation (||) . We have normalized both ex- 
perimental and theoretical data with Ao = 1078 A. The 
agreement between our theory and the data is good at 
B = 4T and 6T, but at B = 0.5T the theory shows less 
linearity at low temperatures compared to the experi- 
mental data, although they agree fairly well at T = 
and high T. 

As we mentioned earlier, a complete calculation should 
include nonlinear corrections jl5,[l6| in addition to the 
nonlocal corrections considered here. However, as ap- 
parent from Fig. at T = the main source of the 
field dependence of A c ff is the nonlocal effect. Since this 
field dependence is closely related to the flattening of the 
curves at low T, it is reasonable to assume that nonlocal 
effects also dominate the T-dependence. The main dis- 
crepancy between the theory and experiment in Fig. |^ is 
the lack of linearity at 0.5T; however it is difficult to en- 
vision how nonlinear corrections could cure this. Given 
the inherent uncertainties in the extraction of A from the 
/iSR data we consider the overall agreement to be rea- 
sonable even without inclusion of the nonlinear effects. 

In summary, we have calculated the field and temper- 
ature dependence of the effective penetration depth A e ff 
from a nonlocal London model of a d-wave supercon- 
ductor. We used a definition of A c ff which permits a 
direct comparison to the /xSR experimental data. Our 
results exhibit a Independence in the X~g curve below 
T* = Ao(£o/cO ~ V^B, quantitatively consistent with 
the experimental data on YBCO |J. This flattening has 
nothing to do with the reduction of the superfluid density 
or opening a true gap at high magnetic fields. Rather, 
it is a consequence of the nonlocal response of a d-wave 
superconductor which modifies the magnetic field distri- 
bution in the vortex lattice as compared to an ordinary 
London model. Thus, as pointed out previously Jl5| , p7t , 
it is essential to make a distinction between the London 



penetration depth Al (as measured e.g. in the microwave 
experiment pTcf] ) and the effective penetration depth A e ff 
deduced from the magnetic field profile in a /xSR experi- 
ment. In the present model A G ff flattens at low tempera- 
tures and finite field while Xl remains linear in T. 
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